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Abstract. A magneto-orbital coupling mechanism is proposed to account for the 
weak energy gap at the Fermi energy in the (tt, 0) ordered SDW state of a realistic 
three band model for iron pnictides involving dxz, dyz, and dxy Fe orbitals. The orbital 
mixing terms between the dxy and dxz/dyz orbitals, which are important in reproducing 
the orbital composition of the elliptical electron pockets at (±7r,0) and (0, ±7r), are 
shown to play a key role in the energy gap formation in the SDW state. 
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The iron pnictides exhibit a rich phase diagram inEi including magnetic, structural 
(tetragonal-to-orthorhombic), and superconducting phase transitions [3]. The magnetic 
state exhibits (tt, 0) ordering of Fe moments in the a-b plane, with a concomitant 
structural distortion a^b possibly correlated with the ferro orbital order dxz>dyz as seen 
in ARPES studies |1]. Inelastic neutron scattering experiments mu have shown well 
dehned spin wave excitations with energy scale ~ 200 meV, persisting even above the 
Neel temperature |^. Evidently, short range antiferromagnetic (AF) and ferromagnetic 
(F) order remain in the a and b directions respectively even above the disordering 
temperature for long-range magnetic order, which may account for the narrow nematic 
phase 019] above the Neel temperature where the ferro orbital order [I] and structural 
distortion survive, as well as the temperature dependence of the measured anisotropies 
in a and b directions of magnetic excitations and resistivity laciEn]. 

Evidently, the complex multi-orbital character of the underlying microscopic 
description of iron pnictides plays an important role in understanding their macroscopic 
behavior as seen from the interplay of magnetic, orbital, structural, and transport 
properties. Another important macroscopic behavior of iron pnictides is associated 
with the energy gap in the (tt, 0) ordered magnetic state, as seen in optical conductivity, 
transport and scanning tunneling microscopy studies; investigating the role of the multi¬ 
orbital character on the magnetic-state electronic structure is therefore of interest. 

For the non-magnetic state, intensive investigations of the electronic structure of 
these materials by hrst-principle calculations HB HU ca El na and angle resolved 
photoemission spectroscopy (ARPES) experiments |T6l El 111 IHl E] have revealed that 
the Fermi surface consists of two nearly circular hole pockets around the center and 
elliptical electron pockets around the corner of the BZ which are primarily contributed 
by dxz, dyz, and dxy Fe 3d orbitals. 

A variety of experimental techniques have confirmed the opening of a small energy 
gap in the excitation spectrum associated with the onset of the (tt, 0) spin-densitywave 
(SDW) order in iron pnictides [20|. The Fermi surface (FS) undergoes a complex 
multi-orbital reconstruction through the non-magnetic to antiferromagnetic transition. 
ARPES experiments have revealed that the electronic states are strongly modihed across 
the SDW transition [211 ^2i l23]. Bands get folded and hybridized, resulting in band 
splitting and opening of SDW gaps. Consequently, the Fermi surface breaks up into 
small droplet-like structures, supporting only a fraction of the original FS areas. These 
results are consistent with Quantum Oscillation measurements [2ll [25] showing the 
reconstruction and drastic reduction of Fermi surface area, which can be as small as 
~ 1% of the original BZ area for SrFe 2 As 2 . 

In addition, investigations of optical and transport properties by optical 
spectroscopy [261 E] and broadband spectroscopic ellipsometry [28] measurements on 
single-crystalline iron pnictides have found an energy gap for electronic excitations in the 
SDW ordered state. The carrier density is reduced substantially in the ordered state with 
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consequent decrease in optical conductivity. Also, as inferred from the results of infrared 
studies, short-range AF correlations have been suggested as being responsible for the 
pseudogap in iron superconductors [291130]. Moreover, scanning tunneling microscopy 
(STM) studies have also seen an energy gap in the density of states near the Fermi 
energy [311132]. 

Close proximity of the superconducting phase and magnetically ordered parent 
phase in iron pnictides have naturally led to comparisons with cuprates [29]. However, 
the magnitude of the experimentally measured energy gap at the Fermi energy (up to 
~ 150 meV for CaFe 2 As 2 ) is signihcantly smaller than the usual SDW gap 2A = mil 
associated with AF ordering with staggered magnetization m and local Coulomb 
interaction U, as in the Hubbard model representation for cuprates. Moreover, these 
SDW gaps open at energies typically away from the Fermi energy in realistic multi¬ 
band tight binding models which yield the correct FS structure in agreement with DFT 
calculations and ARPES studies in the non-magnetic state. 

In this paper, a magneto-orbital coupling mechanism is proposed for a realistic 
three-band model at half hlling for iron pnictides involving d^z, dy^, and d^y orbitals of 
Fe, which is shown to account for the weak non-conventional energy gap at the Fermi 
energy due to a composite effect of the orbital mixing (hybridization) between the d^y 
and dxz/dyz orbitals and the (tt, 0) SDW magnetic order. For realistic three-band-model 
parameters, the energy gap is found to be in quantitative agreement with experiments. 
The gapped SDW state in the realistic three band model at half hlling therefore provides, 
in analogy with cuprates, a suitable reference state for theoretically investigating effects 
of electron/hole doping in iron pnictides. 

2. Three-orbital model and (vr, 0) SDW state 

We consider a minimal three-orbital model [33] involving d^z-, dyz and d^y Fe ?)d orbitals. 
The tight binding Hamiltonian in the plane-wave basis is dehned as: 



( 1 ) 


where 


= — 2ti cos kx — 2^2 cos ky — At^, cos k^ cos ky 
T '22 _ _ 2 -j-^ pQg _ 2 t^ cos ky — 4^3 cos k^ cos ky 
2^33 _ _ 2 ^ 5 (cos kx -l- cos ky) — Ate cos kx cos ky -\- £diff 
= — 4^4 sin kx sin ky 
= — 2^7 sin kx — Aitg sin kx cos ky 
= —2itj sin ky — Ait^ sin ky cos kx 


( 2 ) 


are the tight-binding matrix elements in the unfolded BZ (—vr < kx,ky < vr). Here, 
ti and t 2 are the intra-orbital hoppings for xz {yz) along x {y) and y [x) directions. 


respectively, t^ and ^4 are the intra and inter-orbital hoppings along diagonal direction 
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Figure 1. (a) Fermi surface in the unfolded BZ for the three-orbital model with 


hopping parameters as given in Table [TJ The main orbital contributions are shown 
as: dxz (red), dy^ (green), and d^y (blue). |(b)| The spin wave spectral function in the 
(tt, 0) SDW state for the three-orbital model at half filling. 


Table 1. Values of the hopping parameters in the three-orbital model (in eV) 


h 

^2 


h 

tb 

tb 


ts 

0.1 

0.32 

-0.29 

-0.06 

-0.3 

-0.16 

-0.15 

-0.02 


for xz and yz, t^ and t^ are intra-orbital NN and NNN hoppings for xy, while tj and 
fg the NN and NNN hybridization between xy and xzjyz. Finally, £difr is the energy 
difference between the xy and degenerate xzjyz orbitals. 

The Fermi snrface for the three-orbital Hamiltonian ([1]) for the valnes of hopping 
parameters given in Table 1 is shown in Fignre 1(a). Here, the Fermi energy is kept 
at —0.06 eV corresponding to half hlling. There are two circular hole pockets around 
the center and elliptical electron pockets around (ivr, 0) and (0, ±7r) in the unfolded 
BZ. The two hole pockets involve primarily the xz and yz orbitals, while the electron 
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pockets centered at (ivr, 0) [(0, ±7r)] arise mainly from the hybridization of the xy and 
yz [xz\ orbitals. All of these features are in good agreement with results from DFT 
calculations and ARPES experiments. 

We now consider the (vr, 0) ordered magnetic (SDW) state of this model. The 
various electron-electron interaction terms included are: 

j iJ,<y ^i<u 

HI = U 'y ^ A y ^ 2t/ y ^ ■ Si^j/ 

IJ.<y 

A A H.C.), (3) 

where refer to the local spin (charge) density operators for orbital p. The 

hrst and second terms are the intra-orbital and inter-orbital Coulomb interactions 
respectively, the third term is the Hund’s coupling and the fourth term the “pair¬ 
hopping” term. 

Extending the two-sublattice basis approach for the SDW state in a single-band 
model [21] to a composite three-orbital, two-sublattice basis, the Hartree-Fock (HF) level 
Hamiltonian matrix in this composite basis {Axz Ayz Axy Bxz Byz Bxy) is obtained 
as: 


^hF(k) = 


— crA 


•XZ 1 

0 

0 

4 " + 4 

^k 

4 " + 4 ’ 

0 

-a Ay, + 4^ 

^k 

^k 

4 “+ 4 

4 " 

0 

-e'^y 

-\- £difT 

7x 8,1 

-Ek - ^k 

^k 

-5x _i_ ^6 

^k 

' a4 

^k 

4"+ 4' 

aA^, + 4 ^ 

0 

0 


4 " + 4 


0 

(jAy, + 

^k 

2: ^8.1 


c-Sx 1 -6 

0 

-e'^y 

(TAxy - 




for spin a, where 



- 2ti 

cos kx 

Ay - 

^k “ 

-2ti 

COS ky 


,.2x 

— 2^2 

cos kx 

- 

— 2^2 

cos ky 


^k — 

— 2^5 

cos kx 

Ay — 

^k “ 

— 2^5 

cos ky 


- 

^k — 

— 4f3 

cos kx 

cos ky 

^k = 

—Ati sin kx 

sin ky 

F® - 
^k — 

— Ate 

cos kx 

cos ky 




p.7x _ 

^k — 

— 2it' 

7 sin kx 


—2it 

7 sin ky 


4’^ - 

— Aiti 

i sin kx 

. cos ky 

^k ■ 

= —Aits cos 

kx sin ky 


(5) 


are the band energies corresponding to the hopping terms along different directions, and 
the self-consistent exchange helds are dehned as 2A^ = Um^ J terms of 

sublattice magnetization for orbital y. 
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The calculated spin wave spectral function [33] in the SDW state for the three- 
orbital model is shown in Fig. 1(b). Evidently, spin wave excitations are highly 
dispersive, and do not decay into the particle-hole continuum. The energy scale of 
spin excitations is ~ 200 meV with a well-dehned maximum at the ferromagnetic 
zone boundary IQ = (7r,7r)]. These features of spin wave excitations are in excellent 
agreement with results from inelastic neutron scattering measurements, conhrming that 
our three-orbital model is realistic. 

3. Electronic structure in the SDW state 

Figure 2 shows the the evolution of the orbital resolved density of states (DOS) near 
the Fermi energy in the SDW state with increasing interaction strength U. The 
corresponding band structures along symmetry directions of the BZ are shown in Fig. 
3. Figures 2(a) and 3(a) correspond to the PM state for which the model yields correct 
FS structure as shown in Fig. 1(a). 

As seen in Fig. 2, a large gap opens for the yz band in the SDW state. This 
band gap is a conventional SDW gap, which opens at wave vector fca, = 7r/2 in the AF 
direction, as seen from the band dispersion plot shown in Fig. 3. Similar SDW gaps 
open up for the xz and xy bands at = vr/2. However, these band gaps appear at 
energies far away from the Fermi energy. Besides these three conventional SDW band 
gaps, there is a weak band gap involving the xy and xz orbitals which opens up near 
the Fermi energy, as seen in Fig. 3, which is found to be robust with respect to small 
variations in the hopping parameters. 

It is important to note that while partial density of states for xz and yz orbitals 
are symmetric in the PM state [Fig. 2(a)], they become anisotropic in the SDW state 
[2(b)]. The DOS for the yz orbital has a peak further away from Fermi energy {Ep) 
than that for xz and has a relatively smaller contribution at Ep. This is in excellent 
agreement with results of DMFT studies [35] and provides a good understanding of the 
ferro orbital ordering between the two orbitals [36ll3Z|- Moreover, the band dispersion, 
gap formation, as well as the largest gap for yz [Fig. 2(c) and 2(d)] agree well with the 
self-consistent calculations for multi-orbital models [3H1139] . 

Physically, this weak energy gap near the Fermi energy involving the xy and xz 
orbitals is the most important. One of the key contributing factors for this energy gap 
formation is the orbital mixing between the xz and xy orbitals due to the hopping 
terms ty and ts- Only tj will be retained in the discussion below as tg is negligible in 
comparison (see Table 1). Fig. 4 shows the SDW state band structure with and without 
ty, conhrming the role of the orbital mixing term in the gap formation. It is important 
to note that the opening of this gap for the xz and xy orbitals due to the orbital mixing 
term ty critically requires the presence of the SDW magnetic ordering, as is evident from 
the absence of any gap in the non-magnetic state (A = 0), even with the mixing term 
tj included [Fig. 4 (a)]. For simplicity, we have taken equal exchange helds (A = 0.3 
eV, 17 ~ 1 eV) for all three orbitals. 
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Figure 2. Evolution of the orbitally resolved DOS in the (tt, 0) SDW state of the 
three-band model, showing the opening of the energy gap with increasing interaction 
strength U at the Fermi energy corresponding to half filling. 


To summarize, both the SDW order as well as the orbital mixing term tj between 
the xz and xy orbitals are critically important in the formation of the weak energy gap 
near the Fermi energy. Only the tj term (non-magnetic state) or only the SDW order 
{tj = 0) are independently not sufficient to open the gap, both must be present together. 

4. Magneto-orbital coupling 

A magneto-orbital coupling mechanism is proposed below to account for this non- 
conventional weak energy gap near the Fermi energy involving the xz and xy orbitals. 
In the (vr, 0) ordered SDW state with AF ordering in the x direction, the orbital mixing 
term 2i{t-j + 2ts) sin kx connects xz and xy orbitals on opposite sublattices, highlighting 
the coupling between the orbital and magnetic sectors. 

The composite effect of the orbital mixing between the xz and xy orbitals and 
the (tt, 0) SDW order on opening the weak band gap near the Fermi energy can be 
clearly seen by going over to a reduced basis involving the relevant xz and xy bands and 
contracting over the sublattice basis to obtain the band mixing terms in the magnetic 
state. 

We focus on the kx direction of the Brillouin zone which is the AF ordering direction. 
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Figure 3. Evolution of the band structure in the (tt, 0) SDW state of the three-band 
model with increasing interaction strength U. 



Figure 4. Band dispersion in the (tt, 0) SDW state along the AF direction (F —>■ A), 
showing that the energy gap opens only when both A and are non-zero. The Fermi 
energy corresponds to half filling, for which the FS in the PM state of the three-band 
model is shown in Fig. [I] 
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and set ky = 0. As the yz orbital completely decouples in this case, we start with the 
[4 X 4] part of the Hamiltonian matrix given in Eq. (5) involving the {Axz Bxz Axy 
Bxy) part of the orbital-sublattice basis for the xz and xy orbitals: 



where tk,mix = —‘2.i{ti + 2ts)smkx is the mixing term between xz and xy orbitals on 
opposite sublattices, and the bare band energy terms are e\^^xz = (—2ti — dta) cos kx and 

^k,a;3/ ( 2^5 4^0) COs/Cj,. 

In the absence of the above mixing term tk.mix, the xz and xy orbital sectors are 
decoupled and in the (tt, 0) ordered SDW state, the SDW band energies for the upper 
(-I-) and lower (-) bands are given by: 

^t,xy = “2^5 + Sdiff ± \J + ^"i^xy (V 

The exchange helds Axz and Axy open band gaps at kx = n/2 where £k,x^ = £'k,xy = 0. 
However, these gaps open at energies away from the Fermi energy, as discussed earlier 
[see Fig. 4], For kx close to zero, which is far away from where the band gaps opens, 
the band energies are nearly unchanged for small A [see Fig. 4(c)], and are given by: 

^\l,xz ~ ~2t2 i £'k,xz 

2^5 T A ^\i,xy (8) 

Out of these four bands, we focus on the two relevant SDW bands (upper 

xz band) and E^^y (lower xy band) which he near the Fermi energy. We now consider 
an effective [2 x 2] Hamiltonian matrix in the reduced basis consisting of these two 
eigenstates |k+^) and The band mixing term in this reduced basis is obtained 

by evaluating the matrix elements using the corresponding eigenvectors in the two- 
sublattice basis: 
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is a magneto-orbital coupling term involving the orbital mixing term and the SDW state 
matrix element of representing the opposite sublattice operator. Here ttk and /3k are 
the A and B sublattice components of the eigenvectors for the two orbitals: 



We note here that the majority/minority amplitudes on the two sublattices get 
exchanged for the xz upper band and xy lower band. 

In the absence of any SDW order {A^z = A^y = 0), the magneto-orbital coupling 
term = 0 as the A and B sublattice amplitudes are identical in magnitude 
(Ittkl = |/^k| = 1/\/2), so the negative sign in Eq. (12) yields an exact cancellation. 
However, for hnite and A^y, the sublattice asymmetry in the magnetic state results 
in 7 ^ 0, which leads to strong mixing at k-points where the two band energies 
and are nearly degenerate, resulting in the opening of the non-conventional weak 
energy gap near the Fermi energy. For ~ ^k,xy — ^k — A, we 

obtain: 


rmo 


[ 2 i(t 7 -|- 2t^) sin 



(13) 


explicitly showing the composite dependence of the magneto-orbital coupling term on 
both the orbital mixing and SDW order, and thus accounting for the absence of the band 
gap near the Fermi energy if either of these two terms vanishes, as seen in Fig. 4. 


5. Conclusions 

A magneto-orbital coupling mechanism was proposed and shown to account for the 
opening of a non-conventional weak energy gap near the Fermi energy in the (tt, 0) 
SDW state of a realistic three band model for iron pnictides at half hlling. The coupling 
term involves a composite dependence on both the orbital mixing (hybridization) terms 
between the xy and xzjyz orbitals as well as the SDW order parameter, thus accounting 
for the absence of the energy gap if either of these two terms is absent. As the orbital 
mixing terms also play an important role in the mixed composition {xy and xzjyz') 
of the elliptical electron pockets in the non-magnetic-state FS structure, the magneto¬ 
orbital coupling provides a subtle link between the non-magnetic-state FS features and 
the SDW state weak energy gap. Furthermore, the non-conventional nature of the weak 
energy gap for quasiparticle excitations with orbital mixing provides a sharp contrast 
from the conventional strong SDW gap in cuprate antiferromagnets. 
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